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ABSTRACT 

Context. Solar flares strongly affect the Sun's atmosphere as well as the Earth's environment. Quantifying the maximum 
possible energy of solar flares of the present-day Sun, if any, is thus a key question in heliophysics. 
Aims. The largest solar flares observed over the past few decades have reached energies of a few times 10 32 ergs, possibly 
up to 10 33 ergs. Flares in active Sun-like stars reach up to about 10 36 ergs. In the absence of direct observations of solar 
flares within this range, complementary methods of investigation are needed to assess the probability of solar flares 
beyond those in the observational record. 

Methods. Using historical reports for sunspot and solar active region properties in the photosphere, we scaled to ob- 
served solar values a realistic dimensionless 3D MHD simulation for eruptive flares, which originate from a highly 
sheared bipole. This enabled us to calculate the magnetic fluxes and flare energies in the model in a wide paramater 
space. 

Results. Firstly, commonly observed solar conditions lead to modeled magnetic fluxes and flare energies that are com- 
parable to those estimated from observations. Secondly, we evaluate from observations that 30% of the area of sunspot 
groups are typically involved in flares. This is related to the strong fragmentation of these groups, which naturally 
results from sub-photospheric convection. When the model is scaled to 30% of the area of the largest sunspot group 
ever reported, with its peak magnetic field being set to the strongest value ever measured in a sunspot, it produces a 
flare with a maximum energy of ~ 6 x 10 33 ergs. 

Conclusions. The results of the model suggest that the Sun is able to produce flares up to about six times as energetic 
in total solar irradiance (TSI) fluence as the strongest directly-observed flare from Nov 4, 2003. Sunspot groups larger 
than historically reported would yield superflares for spot pairs that would exceed tens of degrees in extent. We thus 
conjecture that superflare-productive Sun-like stars should have a much stronger dynamo than in the Sun. 

Key words. Magnetohydrodynamics (MHD) - Sun: flares - Solar-terrestrial relations - Stars: flares 



1. Introduction 

Solar flares result from the abrupt release of free magnetic 
energy, that has previously been stored in the coron al mag- 
netic field b y flux emergence and surface motions (Forbes 
et al. 2006). Most of the st rongest flares are eruptive (as 
reviewed by Schrijver 2009). For the latter, the standard 
model attributes the flare energy release to magnetic recon- 
nection that occurs in the wake of coronal mass ejections 



drivers for space weather, and a re responsible for various 
envir onmental hazards at Earth ( Schwenn|2006 Pulkkinen 



2007) 



(IShibata et al.|1995HLin fc Forbes||2000| |Moore et al.|2001 
Priest & Forbes 2002). 



'riest fc Forbes||2002 ) 

Several flare-related phenomena impact the solar at- 
mosphere i tself. To be specific, there are photospheric 
sunquakes (Zharkov et al. 2011), chromospheric ribbons 
dSchmieder et aT. 1987[ ), coronal lo op restructuration 
(Warren et al. 201 ip and oscillation (Nakariakov et al 



For all these reasons, it would be desireable to know 
whether or not there is a maximum for solar flare energies, 
and if so, what its value is. 

On the one hand, detailed analyses of modern data from 
the past half-century imply that solar flare energies range 
from 10 28 to 10 33 ergs , with a power-law d istribution that 
drops above 10 32 ergs ( |Schrijver et al.|2012[ ). The maximum 
value there corresponds to an estimate for the strongest 
directly-observed flare from Nov 4, 2003. Saturated soft X- 
ray observations showed that this flare was above the X28 
class, and model interpretations of radio observations of 



1999D, larg e-scale coronal propagation fronts QDelannee 
et al.|2008 ), and driving of sympathetic eruptions (Schrijver 
fc Title||2011 ). In additi on to solar effects, flare-related ir- 
radiance enhance ments fWoods et al.|2004 ), solar energetic 



Earth ' s iono sphere suggested that it was X40 (Brodrick 
et al.||2005|). Due to the limited range in time of these 



particles (SEPs, j Masson et al 
ejections (CMEs, Vourlidas et a" 



observations, it is unclear whether or not the Sun has 
been -or will be- able to produce more energetic events. 
For example, the energy content of the first-ever observed 



2009) and coronal mass 



solar flare on Sept 1, 1859 ( Carrington] |1859[ |Hod gson 
1859[) has been thoroughly debated QMcOracken et al.J2001 



2TjTo"1 ) constitute major |Tsurutani et al.|2003| [Oliver fc S valgaard |2004| | Wolff et al.| 
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2012|. On the other hand, precise measurements on unre- 
solved active Sun-like stars have revealed the existence of 
so-call ed superflares, even in slowly-rotating and isolated 



it erupted, the flux rope and a surr ounding sigmoid were 



stars (|Schaefer et al 
have 



2000 Maehara et al. 2012) 



Their 

1Q 33 



energies nave been estimated to be between a tew 
ergs to more than 10 36 ergs. Unfortunately, it is still un- 
clear whether or not the Sun can produce such superflares 
among other reasons because of the lack of relia ble informa- 



tion on the starspot pr operties of such stars (Berdyugina 



2005 Strassmeier 2009) 



so as to estimate flare energies, a method complemen- 
tary to observing solar and stellar flares is to use solar 
flare models, and to constrain the parameters using obser- 
vational properties of active regions, rather than those of 
the flares themselves. In the present paper, we perform such 
an analysis. Since analytical approaches are typically over- 
simplified for such a purpose, numerical models are likely 
to be required. Moreover, incorporating observational con- 
straints not only precludes the use of 2D models, but also 
restrict the choice to models that have already proven to 
match various solar observations to some acceptable degree. 

We use a zero-/3 MHD simulation of an eruptive flare 
( Aulanier et al.|2010 2012 ) that extends the standard flare 
model in 3D. Dedicated analyses of the simulation, as re- 
called hereafter, have shown that this model successfully re- 
produced the time-evolution and morphological properties 
of active region magnetic fields after their early emergence 
stage, of coronal sigmoids from their birth to their erup- 
tion, of spreading chromospheric ribbons and sheared flare 
loops, of tear-drop shaped CMEs, and of large-scale coronal 
propagation fronts. We scaled the model to solar observed 
values as follows: we incorporate observational constraints 
known from previously reported statistical studies regard- 
ing the magnetic flux of active regions, as well as the area 
and magnetic field strength of sunspot groups. This method 
allows one to identify the maximum flare energy for real- 
istic but extreme solar conditions, and to predict the size 
of giant starspot pairs that are required to produce super- 
flares. 



progressively formed in the corona ( Aulanier et aL| 2010 
|Savcheva et al.|2012 ), above a slowly shearing and diffusing 



photospheric bipolar magnetic field. This pre-eruptive evo- 
lution was similar to that applied in past symmetric models 



van Ballegooijen & Martens 1989 Amari et al.|2003 l, and 
they matched observations and simulations for active re- 
gions during their late flu x emergence stage and their sub- 



sequent decay phase (e.g. [va n Dricl-Gesztelyi et al. 
Ax iioiiib. et al." 21 )Q4j preen et al.||201ip . 



2003 



The magnetic field geometry of the modeled eruptive 
flare is shown in Fig. [l] The left panel clearly shows the 
asymmetry of the model. A 27% flux imbalance in the pho- 
tosphere, in favor of the positive polarity, manifests itself 
as open magnetic field lines rooted in the positive polarity, 
at the side of the eruption. This asymmetry was set in the 
model so as to reproduce typical solar active regions, with 
a stronger (resp. weaker) leading (resp. trailing) polarity. 
In the right panels, the field of view corresponds to the size 
of the magnetic bipole L bl P ole j as used for physical scaling 
hereafter. 

If one assumes a sunspot field of iJ™ ax = 3500 G, then 
the isocontours that cover the widest areas correspond to 
_BP ax /5 = ±700 G. Since this is th e minimum magneti c 



field value for sunspot penumbrae (Solanki et al. 2006) 



those isocontours correspond to the outer edge of the mod- 
eled sunspots. With these settings, the total sunspot area 
in the model is about half of the area of the field of view 
being shown in Fig. [TJ right. So with Bf ax = 3500 G the 
sunspot area is / _1 (X blpole ) 2 , with / ~ 2, while a lower 
value for £? rnax implies a higher value for /. 

During the pre-eruptive energy storage phase, the com- 
bined effects of shearing motions and magnetic field diffu- 
sion in the photosphere eventually resulted in the devel- 
opment of magnetic shear along the polarity inversion line, 
over a length of about £ bl P ole . This long length presumably 
results in the modeled flare energy to be close to its maxi- 
mum possible value, given the distribution o f photospheric 
flux (|Falconer et al.||2008[ |Moore et al.||2012[). 



2. The eruptive flare model 

2.1. Summary of the non-dimensionalized model 

The eruptive flare model was calculated numerically, using 
the observationally drive n high-order scheme m agnetohy- 
drodynamic code (OHM: Aulanier et al. 2005). The cal- 
culation was performed in the pressureless resistive MHD 
approximation, using non-dimensionalized units, in a 251 x 
251 x 231 non-uniform cartesian mesh. Its uniform resistiv- 
ity resulted in a Reynolds number of about R m ~ 10 3 . The 
simulation setting s are thoroughly described in |Aulanier] 
eTaLl pOlcil [20T2I . 



In the model, the flare resulted from magnetic recon- 
nection occuring at a nearly vertical current sheet, grad- 
ually developing in the wake of a coronal mass ejection. 
The reconnection led to the formation of ribbons and flare 



loops ( Aulanier et al.||2012[ ). The CME itself was triggered 
by the ide al loss-of -cquilibrium of a weakly twisted coronal 
flux rope ( Aulanier et al.|20T0), corresponding to the torus 



instab ility ( |Kliem fc Torok 2006 Demoulin & Aulanier 
2010 1. During the eruption, a coronal propagation front 
developed at the edges of t he expanding sheared arcades 
surrounding the flux rope (Schrijver et al. 2011). Before 



2.2. Physical scalings 

The MHD model was calculated in a wide numerical do- 
main of size 20 x 20 x 30, with a magnetic permeability 
|U = 1, using dimensionless values iJ™ ax — 8 in the domi- 
nant polarity, and £ bl P olc = 5. These settings resulted in a 
dimensionless photospheric flux in side the dominant polar- 
ity of <f> = 42 ( Aulanier et al.|2010 ), and a total pre-eruptive 
magnetic energy of _E bl P° Je = 225. 

Throughout the simulation, a magnetic energy of 

^model = l 9 %£ b ipole = 42 wag re l eased . Only 5% of 

this amount was converted into the kinetic energy of the 
CME. These numbers h ave been presented and discussed 
d2012b. The remaining 95% £ modcl of the 



Aulanier et al. 



magnetic energy release can then be attributed to the flare 
energy itself. 

It must be pointed out that the simulation did not cover 
the full duration of the eruption. Indeed, numerical instabil- 
ities eventually prevented us from pursuing it with accept- 
able diffusion coefficients. Nevertheless, the rate of mag- 
netic energy decrease had started to drop before the end 
of the simulation, and the electric currents within the last 
reconnecting field lines where relatively weak. On the one 
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Size of 
photospheric bipole 




Fig. 1. Eruptive flare model, [left:] 
Projection view of randomly plotted 
coronal magnetic field lines. The grayscale 
corresponds to the vertical component 
of the photospheric magnetic field B z . 
[right:] Photospheric bipole viewed from 
above. The pink (resp. cyan) isocontours 
stand for positive (resp. negative) values of 
Bf^/l.l, 2, 3, 4, 5. The yellow isocontour 
shows the polarity inversion line B z = 0. 
[right-top:] The grayscale for B z is the 
same as in the left panel, [right-bottom:] 
The grayscale shows the vertical compo- 
nent of the photospheric electric currents. 
Strong elongated white/black patches 
highlight flare ribbons. The red lines show 
representative post-reconnection flare 
loops, rooted in the flare ribbons. 



hand, this means that the total energy release £; model is ex- 
pected to be slightly higher, but presumably not by much. 
On the other hand, the relatively low R m value of the sim- 
ulation implies that some amount of i^ model should be at- 
tributed to large-scale diffusion, rather than to the flare 
reconnection. 

Because of these numerical concerns, we consider there- 
after that the flare energy in the model was about E = 40, 
but this number should not be taken as being precise. Also, 
within the pressureless MHD framework of the simulation, 
the model cannot address which part of this energy is con- 
verted into heating, and which remaining part results in 
particle acceleration. 

It is straightforward to scale the model numbers given 
above into physical units. In the international system of 
units (SI), jU = 47iT0~ 7 , the total magnetic flux <j) and the 
total flare energy E can then be written as 



E = 




Wb. 



(1) 



(2) 



Rearranging these equations into commonly used solar 
units leads to: 



(rjmax \ f T bipole 

E = 0.5 x 10- feY fi^V 
I 10 3 G \ 50 Mm j 



Mx. 



erg . 



(3) 
(4) 



While the power-law dependences in these equations come 
from the definitions of flux and energy, the numbers them- 
selves directly result from the MHD simulation, and not 
from simple order of magnitude estimates. So Eqs. ^ and 
Q enable us to calculate the model predictions for a wide 
range of photospheric magnetic fields and bipole sizes. The 
results are plotted in Fig. [2} In this figure, the right vertical 
axis is the total sunspot area within the model, being given 



using / = 2. It is expressed in micro so- 
lar hemispheres (hereafter written MSH as in Baumann & 
Solanki 2005, although other notations can be found in the 
literature). Hereafter all calculated energies (resp. fluxes) 
will almost always be given in multiples of 10 32 ergs (resp. 
10 22 Mx), for easier comparison between different values. 



Typical decaying active regions with L bi P° lc = 200 Mm, 
which contain faculae of Bf ax = 100 G, have <j> = 0.8 X 10 22 
Mx and can produce moderate flares of E — 0.3 x 10 32 ergs. 
Also, J-spots with L bi P° le = 40 Mm and 5™ ax = 1500 G 
have a lower magnetic flux (f> — 0.5 x 10 22 Mx, but can 
produce twice stronger flares, with E = 0.6 x 10 32 ergs. 
These energies for typical solar active regions are in good 
agreement with those estimated from the tota l solar irradi- 
ance (TSI) fluencc of several observed flares ( Kretzschmar 
20111. 



Other parameters can result in more or less energetic 
events. For example one can scale the model to the sunspot 
group from which the 2003 Halloween flares originated. 
Firstly, one can ove rplot our Fig. [Tj right , onto the cen- 



ter of the Fig. 2 in Schrijver et al. (20061) and thus find 



an approximated size of the main bipole which is involved 
in the flare, out of the whole sunspot group. This gives a 
bipole size of the order of £ bi P° le ~ 65 Mm. Secondly, ob- 
servational records lead to a peak sunspo t magnetic field 
of Bf ax = 3500 G ( [Livingston et~alp012l) . These scalings 
lead to (j) = 3 x 10 22 Mx and E = 13x 10^ ergs. The mod- 
eled 4> is about one third of the flux of the dominant polarity 
as me asured in the whole active region (Kazachenko et al. 
2010). Comparing this modeled flare energy E with that 
of extreme solar flares that originated from this same ac- 
tive region, we find that it is twice as strong a s that of the 
Oct 28, 2003 X17 flare flSchrij ver et al.]|2012| ) , and about 
the same as that of the Nov 4, 2003 X28-40 flare, as can 



be est imated from Kretzschmar (2011) and Schrijver et al. 
p012l Eq. 1). 
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Flux of the strongest polarity (Mx) Energy of the flare (erg) 




10 100 1000 10000 10 100 1000 10000 

Max of photospheric B (G) in bipole Max of photospheric B (G) in bipole 



Fig. 2. Magnetic flux in the dominant polarity of the bipole, and magnetic energy released during the flare, calculated 
as a function of the maximum magnetic field and the size of the photospheric bipole. The x and + signs correspond to 
extreme solar values. The former is unrealistic and the latter must be very rare (see text for details). 



3. Finding the upper limit on flare energy 

3.1. Excluding unobserved regions in the parameter space 

We indicate in Fig. [2] the minimum and maximum sunspot 
magnetic fields as measured from spectro-polarimetric ob- 
servations since 1957. They are respe ctively 700 G in the 



penumbra, and 3500 G in the umbra (Solanki et al. 2006 



Pevtsovet al.|2011[ ). The latter value is an extreme that has 



rarely been reported in sunspot observations, and it typi- 
cally is observed in asso ciation with intense flaring activity 
( Livingston et al.||2012 ). 

We also indicate the maximum area of sunspot groups, 
including both the umbras an d the penumbras. They were 
measured from 1874 to 1 976 flBaumann fc Solanki||2 005) 
and from 1977 to 2007 QMathaway fc Ghoudhary] |2008 ) . 
These sizes follow a log-normal distribution up 3000 MSH, 
but there are a few larger groups. The largest one was ob- 
served in April 1947, and its area w as about 5400 — 6000 
MSH {Nicholson! 119481 |Taylor||1989b. For illustration, we 



provide in Fig. [3 one image of this sunspot group and one 
of its surrounding faculae and filaments, as observed with 
the Meudon spectroheliograph. Interestingly, this sunspot 
group did not generate strong geomagnetic disturbances. 
This could either be due to a lack of strong enough mag- 
netic shear in the filaments which were located between the 
sunspots, or to the lack of Earth-directed CMEs that could 
have been launched from this region. However, several other 
large sunspot groups, whose areas were at least 3500 MSH, 
did generate major geomagnetic storms. Among those are 
the March 1989 event, which led to the Quebec black- 
out (Taylor 1989), and the December 1128 event, which 
produced aurorae in Asia and which corresponds to th e 
first reported sunspot drawing dWillis fc Stephenson|206T| ). 
Therefore, we conservatively keep 6000 MSH as the maxi- 
mum value. The 1874-2007 dataset does not include the first 
observed flare, in December 1859. Nevertheless, Hodgson 



mate from the drawing of Carrington ( 1859[ ) that its total 
area was smaller than 6000 MSH. 

The point marked by a thick x sign in Fig. [2] is defined 
by the intersection of the 3500 G and the 6000 MSH lines. 
The model states that its magnetic flux is <f> — 27 x 10 22 
Mx. This modeled value is much higher than 8 x 10 22 Mx, 
which correspon ds both to the dominant polarity for the 
Halloween flares ( Kazachenko et aL]|2010 | and to the high- 
est flux measured for single active regions, as observed 
during a samp l e of t ime-periods between 1998 and 2007 
(Parnell et al. 2009). The modeled flux for this largest 



sunspot group is nevertheless consistent with the maximum 
val ue of 20 x 10 22 Mx for an active region, as reported 
by Zhang et al. (2010) in a very extensive survey, rang- 
ing from 1996 to 2008. It remains difficult to estimate the 
highest active region flux which ever occurred. Firstly, no 
magnetic field measurement is available for the April 1947 
sunspot grou p. Secondly, the automatic procedure of |Zhang] 
et al. (2010) can lead several active regions to be grouped 
into an app arent single region, while the method of |ParneII| 
et al. ( 2009 1 in contrast tends to fragment active region into 
several pieces. For reference, we therefore overplotted both 
(f> = 8 x 10 22 and 2 x 10 23 Mx values in Fig. [2 

The flare energy at the point x , where the magnetic 
field and size of sunspot groups take their extreme values, 
is E x — 340 x 10 32 ergs. This could a priori be considered as 
the maximum possible energy of a solar flare. In a ddition, it 



( 1859 ) reported that the size of the sunspot group associ- 
ated with this event was about 96 Mm, and one can esti- 



falls within the range of stellar superflare energies (Maehara 
et al. 2012). Nevertheless, we argue below that this point 
is unrealistic for observed solar conditions. 



3.2. Taking into account the fragmentation of flux 

All large sunspot groups are highly fragmented, and display 
many episodes of flux emergence and dispersal. We argue 
that this fragmentation is the reason why scaling the model 
to the whole area of the largest sunspot group leads to over 
estimate the maximum flare energy. 
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Fig. 3. The largest sunspot group ever reported since the end of the nineteenth century, as observed in April 5, 1947 in 
Ca II Klv (left) and Ha (right) by the Meudon spectroheliograph. 



Firstly, sunspot groups incorpo rate several big sunsp ots 
ranging from a few spots (see e.g. Schrijver et al.||2011[ for 
February 2011) to half a dozen (see e.g. Carrington 1859 



|Schrijver et al.|2006[ for September 1859 and October 2001 



respectively) and up t o more than ten (see e.g. Wang et al. 
1991] |Nicholson|[l94"8l for March 1989 and April 1947 re- 
spectively; see also Fig. [3]). Secondly, these groups typically 
have a magnetic flu x imbalance (e.g. 23% fo r the October 
2003 sunspot group Kazachenko et aL]|2010 1 , because they 
often emerge within older active regions. This naturally cre- 
ates new magnetic connections to distant regions on the 
Sun, in addition to possibly pre-existing ones. Thirdly, the 
magnetic shear tends to be concentrated along some seg- 
ments only of the pol arity inversion lines of a given group 
( Falconer et al.||2008| ). This is also true for the April 1947 
sunspot group, as evidenced by the complex distribution 
of small filaments (see Fig. [3]). This means that a given 
sunspot group is never energized as a whole. These three 
observational properties are actually consistent with the 
solar convection-driven breaking of large sub-photospheric 
flux tubes into a series of smal ler deformed structures, as 
found in numerical simulations ( |Fan et al.|2003 Jouve et al. 
20131. They show that these deformed structures should 



eventually emerge through the photosphere as grouped but 
distinct magnetic bipoles. These different bipoles should 
naturally possess various degrees of magnetic shear, and 
should not be fully magnetically connected to each other in 
the corona. 

So both observational and theoretical arguments sug- 
gest that only a few sunspots from a whole sunspot group 
should be involved in a given flare. Unfortunately, the frac- 
tion of area to be considered, and to be compared with the 
size of the bipole in the model, is difficult to estimate. 

We consider the Oct-Nov 2003 flares, for example. Our 
estimation of L hi P ole ~ 65 Mm, as given above, results in a 
modeled sunspot area of 700 MSH (see Fig.[2|. This is about 
27% of the maximum area measured for the whole sunspot 
group, which peaked at 2600 MSH on Oct 31. Another way 
to estimate this fraction is to measure the ratio between 



the magnetic flux swept by the flare ribb ons, and that of 
the whole ac tive region. Qiu et al. ( 2007| and |Kazachenko 
et al.| ( |2012[ ) reported a ratio of 25%and 31% for the Oct 
28 flare, respectively. The same authors also reported on 
a dozen of other events, for which one can estimate ratios 
ranging between 10% and 30%, on average. 

These considerations lead us to conjecture that 30% at 
most of the ar ea of the l argest obse rved sunspot group, 
as reported by |Nicholson| ( |1948| and |Taylo"r] ( |1989| ), i.e. a 
maximum of 1800 MSH, can be involved in a flare. This is 
more than 2.5 times the area of the bipole involved in the 
Halloween flares. In Fig. [2j we therefore plot another point 
indicated by a thick + sign, located at the intersection of 
the 3500 G and the 1800 MSH lines. In the model, this 
corresponds to £ bl P olc = 105 Mm. The flare energy at this 
point is E + = 56 x 10 32 ergs. Under the assumptions of 
the model, and considering that it probably corresponds 
to the most extreme observed solar conditions, E + should 
correspond to the upper limit on solar flare energy. 



3.3. Numerical concerns 

As for all numerical models, various limitations could play 
a role in changing the estimated maximum flare energy E + . 

We mentioned above that the simulation did not cover 
the duration of the full eruption, because some numerical 
instabilities eventually developed. On the one hand, this 
means that our flare energies are slightly under estimated. 
But on the other hand, the low R m must lead to a weak 
large-scale diffusion. It should not be very strong, however, 
since the characteristic diffusion time at the scale of the 
modeled bipole can be estimated as 150 times the dura- 
tion of the simulation. Still, it ought to take away some 
fraction of the magnetic energy released during the simula- 
tion, so that our flare energies are slightly over estimated. 
Quantifying the relative importance of both effects is un- 
fortunately hard to achieve. 

Moreover, applying different spatial distribution of 
shear during the pre- flare energy storage phase could lead to 
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a different amount of energy release ( Falconer et al.| [2008 ) . 
But in our model, the shearing motions were extended ail 
along the polarity inversion line in the middle of the flux 
concentrations. Therefore we argue that it will be difficult 
for different settings to produce significantly higher flare 
energies. 

Another concern is that our simulation produces a CME 
kinetic energy which is only 5% of the flare energy. But cur- 
rent observational energy estimates im ply that the kinetic 
energy of a CME can be the same as (lEmslie et al.||2005| ) 
and up to three times higher than ( |Emslie et al.|2*012| ) the 
bolometric energy of its associated flare. This strong dis- 
crepancy cannot be attributed to the fact that our simu- 
lation was limited in time. Indeed, other 3D (resp. 2.5D) 
MHD models calculated by independent groups and codes 
predict that no more than 10% (resp. 30%) of the total 
released mag netic energy is converted into the CME ki 
netic energy (|Amari et al. 2003 



2t 



Jacobs et al. 2006 Lynch 



et al. 2008 Reeves et al. |2010[ ). This means that it is un- 
ciear whether the relatively weaker CME kinetic energy in 
our model should be attributed to observational biases, or 
to numerical problems commonly shared by several groups 
and codes. 

In principle, the validity of the model can also be ques- 
tioned because magnetic reconnection is ensured by resis- 
tivity, with a relatively low magnetic Reynolds number R m 
as compared to that of the solar corona. This may lead to 
different reconnection rates from those found in collision- 



less reconnection simulations (see e.g. Aunai et al. 2011). 
The reconnection rate is indeed important for the flare en- 
ergy release in fully three-dimensional simulations of solar 
eruptions. In principle, slower (resp. faster) reconnection 
releases weaker (resp. larger) amounts of magnetic energy 
per unit time. Nevertheless, one might argue that the time- 
integrated energy release, during the whole flare, could be 
not very sensitive to the reconnection rate. However the 
energy content which is available at a given time, within a 
given pair of pre-reconnecting magnetic field lines, strongly 
depends on how much time th ese field lines have h ad to 
stretch ideally (as described in Aulanier et al. 2012), and 
thus by how much their magnetic shear has decreased be- 
fore they reconnect. This explains why the time-evolution 
of the eruption makes the reconnection rate important for 
time-integrated energy release. In our simulation, we mea- 
sure the reconnection rate from the average Mach num- 
ber M of the reconnection inflown. During the eruption, 
it increases in time from M ~ 0.05 to M ~ 0.2 approx- 
imately. These reconnection rates are fortunately compa- 
rable to those obtained for collisionless reconnection. So 
we conjecture that the limited physics inside our modeled 
reconnecting current sheet should not have drastic conse- 
quences for the flare energies. Nevertheless, it should be 
noted that this result probably does not hold for other re- 
sistive MHD simulations that use very different R m . 

We foresee that these numerical concerns are probably 
not extremely sensitive: the orders of magnitudes that we 
find for flare energies are likely to be correct. But it is dif- 
ficult at present to assert that we estimate flare energies 
with a precision better than several tens of percents, or 
even more. Therefore we conservatively round up the upper 
value E + to 6 x 10 33 ergs. In the future, data-driven sim- 
ulations which can explore the parameter space and which 
incorporate more physics will have to be developed to fine- 
tune the present analyses. 



4. Summary and discussion 

So as to estimate the maximum possible energy of solar 
flare, we used a dimensionless numerical 3D MHD simu- 

We 



2012) 



lation for solar eruptions (Aulanier et al. 2010 
had previously shown that this model successfully matches 
the observations of active region magnetic fields, of coro- 
nal sigmoids, of flare ribbons and loops, of CMEs, and of 
large-scale propagation fronts. 

We scaled the model parameters to physical values. 
Typical solar active region parameters resulted in typically 



obser ved magnetic fluxe s (|Parnell et al.||2009[ |Zhang et al. 
and flare energies ( |Kretzschmar|2011||Schrijver et al 



2010 



2012). We then scaled t he model using the l argest measured 



sunspot magnetic field ( Solanki et al.|2006 ), and the area of 
the largest sunspot g roup ever reported, which de veloped 
in March- April 1947 ( |Nicholson|[l948"l |Taylor|i~989"l ) . 

In addition, we took into account that observations show 
that large sunspots groups are always fragmented into sev- 
eral spots, and are never involved in a given flare as a 
whole. This partitioning can presumably be attributed to 
sub-photospheric convective motions. Since those motions 
are always present bec ause of the solar internal structure 



2002), it is difficult to imagine that the 



(Brun & Toomre 

Sun will ever produce a large sunspot group consisting of a 
single pair of giant sunspots. Based on some approximated 
geometrical and reconnected magnetic flux estimations, we 
considered that only 30% the area of a given sunspot group 
can be involved in a flare. 

Keeping in mind the assumptions and limitations of the 
numerical model, these scalings resulted in a maximum flare 
energy of ~ 6 x 10 33 ergs. This is is ten times the energy of 
the O ct 28, 2003 X17 flare, as reported in |Schrijver et al.| 
(2012). In addition, this value is about six times higher 
than the maximum energy in TSI fluence that can be es- 
timated from the SXR fluence of the Nov 4, 2003 X2 8-40 
flare, using the scalin gs given by Kretzschmar (20111 and 



Schrijver et al. (2012). Finally, it lies in the energy range 
of the weak est superflares that were reported by Machara 



et al. ( 2012h for numerous slowly-rotating and isolated Sun- 



iike stars. But it is several orders of magnitude smaller than 
that of strong stellar superflares. 

One could ask what the frequency is at which the 
Sun can produce a maximum flare like this. Observational 
records since 1874 reveal that the ar ea of sunspot groups 



follow a sharp log-normal distribution ( Baumann & Solanki 
120051 iHathaway & Choudharyl 120081). Unfortunately, the 



statistics for sunspot groups larger than 3000 MSH in area 
are too poor to estimate whether or not this distribution is 
valid up to 6000 MSH. In addition, neither do all active re- 
gions or sunspot groups generate flares, nor do they always 
generate them at the maximum energy, as calculated by the 
model. The reason must be that a solar eruption requires 
a strong magnetic sheared polarity inversion line, and cur- 
rent observations show that this does not occur in all so- 
lar active regions ( Falconer et al.|[2008[ ). Consequently it is 
currently difficult to estimate the probability of appearance 
of the strongest flare that w e fo und. We can only refer to 



Bau mann fc Solanki (2005) and Hathaway & Choudhary 
( 2008 ), who reported that the size of sunspot groups fol- 



lo ws a clear log-normal distribution up to 3000 MSH, an d 



to 



Cliver k Svalgaard (2004) and Schrijver et al. (2012), 



who argue that this upper limit on flare energy was never 
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produce, a flare more energetic than this upper limit. We 

thus conjecture that one condition for Sun-like stars to pro- 

10 32 erg duce superflares is to host a dynamo that is much stronger 

■ than that of an aged Sun with a rotation rate exceeding 

10 JJ erq , , 

several days. 

On the one hand, our results suggest that we have not 
experienced the largest possible solar flare. But on the 
other hand, and unless the dynamo theory proves other- 
wise, our results also provide an upper limit for extreme 
10 35 erg space weather conditions, that does not exceed those re- 

lated to past observed flares by much. 
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from the AXA Research Fund. 




Fig. 4. Schematic representation of several modeled 
sunspot pairs on the solar disk, with their corresponding 
modeled flare energies. Note that our estimations state that 
in the real Sun, a given pair will often be embedded in a 
much larger sunspot group, from which only the bipole that 
is shown here will be involved in the Hare. 



reached in any observed solar flare since, even including the 
Carrington event of Sept 1859. 

When the model is scaled to the strongest measured 
sunspot magnetic field, i.e. 3.5 kG, it can be used to calcu- 
late the size of the sunspot pair that is required to gener- 
ate the solar flares of various energies. We plot those in 
Fig. |4j These scalings can also be used to relate stellar 
superflares to starspot sizes. But it should be noted that 
starspot magnetic fields are still difficult to measure reli- 
ably, and that current estim ates put them in the range of 
2 — 5 kG (Berdyugina 2005). With these scalings, a super- 
flare of 10 ergs requires a very large single pair of spots, 
whose extent is 48° in longitude/latitude, at the surface of 
a Sun-like star. While such spots have been observed indi- 
rectly in non-S un-like stars as wel l as in young fast-ro tating 
Sun-like stars (Berdyugina 2005 Strassmeier 2009), they 
have never been reported on the Sun. 



5. Conclusion 

We combined a numerical magnetohydrodynamic model for 
solar eruptions calculated with the OHM code and histori- 
cal sunspot observations starting from the end of the nine- 
teenth century. We concluded that the maximum energy 
of solar flares is about six times that of the strongest-ever 
directly-observed flare of Nov 4, 2003. 

One unaddressed question is whether or not the current 
solar convective dynamo can produce much larger sunspot 
groups, as required to produce even stronger flares accord- 
ing to our results. This seems unlikely, since such giant 
sunspot groups "have not been recorded in four centuries 
of direct scientific observations and in millennia of sunrises 
and sunsets viewable b y anyone around the world" , to quote 
Schrijver et al. (2012 1. It can thus reasonably be assumed 
that, during the most recent few billion years while on the 
main sequence, the Sun never has produced, and never will 
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